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ABSTRACT 

Damage  sustained  by  equipment  used  in  a  moving  carrier 
suddenly  changing     its  velocity  can  sometimes  be  attributed  to  in- 
adequate protective  cushioning.     Proper  design  of  the  protective 
cushioning  can  be  realized  by  reducing  the  equipment  and  components 
to  idealized  mechanical  systems  and  studying  the  responses  as  the 
system  parameters  are  varied.    This  study  considers  the  response  to 
velocity  shock  with  elastic  impact.    A  set  of  design  curves  showing 
maximum  acceleration  of  a  flexible  and  vulnerable  element  of  an 
equipment  and  the  maximum  deflection  of  the  equipment's  cushioning 
are  produced  using  an  analog  computer.    Additionally  the  computer's 
results  are  compared  with  analytical  results  for  a  few  particular 
cases . 

The  writer  wishes  to  express  his  appreciation  for  the  advice 
and  assistance  of  Professor  R.  E.  Newton  in  the  selection  and  per- 
formance of  this  study . 
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I.         Introduction 

One  of  the  causes  of  mechanical  damage  in  the  transportation  of 
packaged  articles  is  inadequate  protective  cushioning.    However,  it  is 
possible  to  predict  what  cushioning  is  required  to  protect  the  flexible 
and  vulnerable  elements  of  a  given  equipment.    This  prediction  is  based 
on  the  supposition  that  the  container  for  the  equipment,  the  protective 
cushioning ,  the  equipment  and  the  elements  of  the  equipment  can  be 
idealized  as  a  system  of  rigid  masses,  massless  linear  springs  and 
linear  viscous  dampers  (dashpots) . 

Figs .  1  and  2  illustrate  the  schematic  representation  of  a  package 
and  its  representative  mechanical  system. 


Fig.  1.    Schematic  representation  of  a  package 
1.        Element  of  packaged  item 
2a .      The  packaged  item  as  a  whole 
2b.       Cushioning 
3.        Outer  container 


Element  of  packaged  item 


Packaged  item 


2.  Cushioning 

Outer  container 


\A\\\\\\\\\\\\\\\\\    Support 
Fig.  2.    Idealized  mechanical  system  representing  a  package. 


In  this  study  we  shall  be  concerned  with  the  responses  of  the 
package  that  is  undergoing  velocity  shock  with  elastic  impact*    Velo- 
city shock  is  that  shock  which  results  from  a  sudden  change  in  the  linear 
velocity  of  the  equipment  or  its  support.    If  the  outer  container  nig  does 
not  rebound  from  the  support  we  call  this  behavior  inelastic  impact.    On 
the  other  hand  if  mass  m^  is  negligibly  small  it  will  separate  from  the 
support  when  the  net  force  applied  by  spring  k«  and  dashpot  c«  becomes 
tensile.    We  call  this  latter  behavior  elastic  impact,  or  perfect  rebound. 
The  idealized  system  we  shall  use  exhibits  perfect  rebound  and  is  shown 
in  Fig.  3. 


[fir  I  sw         qq  :  i  • 

'  ■  ■ 

■  ! 

r.rdii  "  -  Gsbi  eri 


a±     t  =  0   : 

o  o  o 

A,  =  Xz.  -  A0. 

WVWWWWWvWWW 

Fig.  3.     Idealized  mechanical  system  consisting  of  a  massive  packaged 

item  m^  and  a  small  element  m^  .    The  outer  container  is  assumed 
massless . 

We  now  have  a  system  with  two  degrees  of  freedom  and ,  to  make 
this  system  tractable,  we  will  consider  it  to  be  linear.    The  free  body 
diagrams  of  the  two  masses  are  shown  in  Fig.  4  for  u  =  u  =  o. 


,ki(X,-Xa) 


k,(x,-xji 

m* 

.   K«.  Xst 

c,(x,-xJi 

c*X* 

Fig.  4.     Free  body  diagrams  for  m,  and  m    of  the  idealized  mechanical 
system. 


The  differential  equations  of  motion  for  m^  and  m„  are: 


YY\ 


.  x,  =  -k,(x,-  xj  -   c,(  x,  -  xj. 


(1) 


>Y\a.X*=  k,(XrXx)+  c,(X,-X^-KxX,-C^XA.    (2) 


It  can  be  seen  that  there  are  six  parameters:    m,  ,  m    ,  k    ,  k„ ,  c",  and  C2 


. 


' 


It  is  convenient  to  define 


CO 
I 


— i 
=    Jii.  CJ„  = 


L-  ,  a  =  -£i      >  ft  =  _£^ 


where  the  new  parameters  have  the  meanings: 

U) ,  ,  the  natural  undamped  frequency  of  vibration  of  mass  m.  on  spring  k.. 
(J  « /  the  natural  undamped  frequency  of  vibration  of  mass  i^  on  spring  k . 
Q>  l ,  the  fraction  of  critical  damping  for  the  subsystem  consisting  of 

ml '  ^1 '  cl  • 

P2#  the  fraction  of  critical  damping  for  the  subsystem  nu ,  k0),  c^. 

Because  both  the  length  and  time  units  may  be  chosen  arbitrarily, 
it  is  possible  to  ignore  absolute  values  of  masses  ihu  ,  m^ ,  and  frequen- 
cies CO  i ,  ^2*    ^e  use  instead  the  mass  ratio    IX  -  mj/m^  and  the  fre~ 
quency  ratio      0J{  /(a)*  •    The  differential  equations  of  motion  become: 
X,    =    ^cot(Xi-Xx)    -^,w,  C*.  -xa).  (4) 

X^     =    -6JA    X^    -  ^a  Wx  X4    -/CX,  .  (5) 

For  perfect  rebound  the  first  two  terms  on  the  right  side  of  Eq.  (S)  do  not 
apply  after  loss  of  contact  with  the  support,    Eq.  (5)  then  becomes: 

5L  =  -ju,  X,  C5-5 

From  a  practical  standpoint  ^  is  normally  small  and  y«>  is 
usually  small.    If  we  fix  fi}  at  a  small  value,  say  0,01,  and  use  the 
limiting  form  assumed  by  Eq.  (5)  as  J/L  approaches  zero,  we  can  reduce 
the  number  of  system  parameters  to  two:        42>  a         and      OJt  /  (O^    . 
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The  behavio      of  the  system  thus  idealized  is  studied  in  Section  II  and 
Appendix  B . 

A  second  two-parameter  idealization  results  from  neglecting 
damping  [pt-B-OJ  and  confining  attention  to  the  effects  of  mass  ratio   M. 
and  frequency  ratio    W,/(aJx  .    This  system  is  studied  in  Section  II  and 
Appendix  B . 

With  just  two  parameters,  analytic  solutions  by  known  methods 
are  feasible.    However,  such  solutions  are  so  laborious  that  it  is  im- 
practical to  produce  a  set  of  design  curves  covering  a  range  of  values 
of  these  parameters.    Accordingly,  a  set  of  design  curves  was  obtained 
using  an  analog  computer. 

0  0 

The  design  results  include  the  maximum  acceleration  X,  ma* 
of  the  element  W  of  the- packaged  item  and  the  maximum  deflection 
of  the  cushioning.    The  maximum  acceleration  of    YV\     is  of  interest 

CO 

because  it  is  a  direct  measure  of  the  maximum  force       \Y[{  X|rAAX  acting 
on  the  structure  of  the  element.    The  maximum  deflection    Xa-wxax      is 
of  interest  because  it  determines  the  required  physical  dimensions  of 
the  cushioning. 


II  .      Experimental  Procedure 

The  free  body  diagrams  for   w,    and    va,    of  the  idealized  mechanical 
system  shown  in  Fig.  4  of  the  Introduction  are  repeated  below. 


k,lx.-xJ 


(i-xj 


ki(x.-xj 

.    Kfl-       ^ 

>*•*. 

dUrXJ. 

c*  *„_ 

at      t  =  0  : 
X,  =  X.  =  0    3 

x,.=  x,=o. 


Fig.  4.    Repeated. 

The  differential  equations  of  motion  are  Eqs .  (4)  and  (5)  of  the 
Introduction . 


By  setting: 


5,  =    X»  "  *a.; 


(5) 


(6) 


; 


(7) 


Eq.  (4)  becomes: 

Eqs.  (5)  and  (7)  are  set  up  on  the  analog  computer  as  indicated  in 
Fig.   5.    Equipment  used  is  listed  in  Appendix  E. 

The  computer  was  scaled  on  a  real  time  basis.    Coefficient 
potentiometers  were  used  to  correct  for  the  differences  between  nominal 


oX 


c    o 
-   u 


i° 


•J 


^ 


X 


and  actual  values  of  resistors  and  capacitors.    The  Donner  Comparator 
Bridge  was  used  for  precision  measurement  of  the  potentiometer  settings. 
A  diode  limiter,  type  6AL5  ,  was  connected  across  the  feedback  ol   3 
summing  amplifier  to  cut  off  the  output  voltage  when  the  voltage  changed 
sign  from  minus  to  plus.    The  time  of  cutoff  is  the  time  of  rebound. 

In  making  the  "runs"  the  initial  condition  voltage    Cjl)^a0      was 
connected  across  the  #3  operational  amplifier.    For  the  ""runs"  with 
damping  and  zero  mass  ratio,  the    Zpt    coefficient  potentiometer  was  set 
at  0.02,  the    2-p%.    coefficient  potentiometer  was  set  as  required  (0.0, 
0.1,  0.2,  0.4,  0.6  or  0.8)  and  the    u    coefficient  potentiometer  was  set 
at  zero.    For  the  "runs"  with  no  damping  and  significant  mass  ratio,  the 
?.£/    and    2.$^  coefficient  potentiometers  were  set  at  zero  and  the    M- 
coefficient  potentiometer  was  set  as  required  (0.01,  0.1,  0.2,  0.4  or 
1.0).    In  all  "runs"  the  Oz  coefficient  potentiometers  were  kept  at  a 
constant  value  of  0.1  rad./sec.  and  the   COl    coefficient  potentiometers 
were  varied  to  provide  the  desired  frequency  ratio  (    CJ,/  60^  ) . 


III.      Analog  Results  for  Velocity  Shock  with  Perfect.  Rebound 

The  analog  results  are  presented  in  the  form  of  design  curves  , 
Figs.  6  to  13.    The  curves  are  plots  of: 

(a)  Xi  ,w;u/  Xt  COv  '  with  frequency  ratio   (*),/&%     as  abscissa. 
Xt  max  /^otJj                is  dimensionless  and  is  called  the  "shock 

transmissibility"  .    It  is  the  ratio  of  maximum  acceleration      X%  ma^ 
to  the  maximum  acceleration     X0  CO,        which  would  result  with  no 
cushioning. 

(b)  Xt  max/  Ac  (j^x  >  with  frequency  ratio  k)JW*   as  abscissa. 
Xi  ma*  /  %o  6J*               is  called  an  "amplification  factor".    It  is  the 


ratio  of  the  maximum  acceleration  of  m,  to  the  maximum  acceleration 
of  m^  with  no  damping  in  the  cushioning. 

(°)    Xa  max/(  Ao/6Ja)       ,  with  frequency  ratio      Ljju>x       as  abscissa. 
V  //  Z    i^     }         is  called  a  "deflection  factor".    It  is  the  ratio 

of  the  maximum  deflection  of  the  cushioning  to  the  rrkaximum  deflection 
of  the  cushioning  with  no  damping  and  negligible  mass  of  the  element  m.  . 
(d)    Xa.  x*\i.y(%0 / to*.)       with  damping  fraction  fi^   as  abscissa.    For  mass 
ratio     >tc  =  °  ,       Xa  m3x/(^0/^z)      is  independent  of  the  frequency 
ratio  but  is  dependent  on    fex      which  can  be  seen  from  Eq»  (22). 

These  curves  can  be  utilized  in  the  selection  of  adequate  cushion- 
ing for  a  flexible  and  vulnerable  element  of  equipment  undergoing  a 
velocity  shock. 


Fig.  7  .  Shock  Transmissibility  for  linear  undamped  cushioning 

and  significant  mass  ratios  with  perfect  rebound 
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Fig.  8  .  Shock  Transmissibility  for  linear  damped  cushioning 

with  perfect  rebound.    Analog  values  are  shown. 
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Fig.  9.  Shock  Transmissibility  for  linear  undampe  I  cushionin 

and  significant  mass  ratios  with  perfec*:  rel  o  ir 
Analog  values  are  shown. 
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Fig.  11.        Amplification  Factors  for  undamped  cushioning  and       ■  :..;  <   tnt 
mass  ratios  with  perfect,  rebound 


16 


Fig 


Def  lectionFactors  for  undamped  cushioning    ,nd  sign 

mass  ratios  with  perfect  reboun  i 


IV.       Comparison  of  Analog  and  Analytical  Resull 

The  analytical  results  of  Appendix  C  are  comp  wit!      r  ilo  | 

results  in  Figs.   14  to  18. 

Figures  16  and  17  show  the  response  cur   e     foi  two  ran 
made  with  the  same  values  for  the  parameters.    There  is    •    lifference  oi 
approximately  two  per  cent  in  the  maximum  magnitu  le     oi  tl     two  runs 

In  all  cases  the  analog  plot  lags  the  analytical  point  plot0    1  \  e 
lag  varies  from  plot  to  plot  and  ranges  from  one  to  three  seconds  In 
seconds.    It  is  difficult  to  determine  a  precise  cor  ,  'ween  r    i  | 

nitudes  of  the  analog  and  analytical  results.    Thi  •        r  .    Is  due  to 

the  time  lag  and  to  the  uncertainty  of  the  exact  analytic  al  maximum  mag- 
nitude.   A  comparison  of  the  maximum  magnitudes  of  the  analog  results 
with  the  largest  analytical  values  suggests  that  the  erroi  r    within  tl  e 
range  of  plus  or  minus  four  per  cento 

A  second  method  of  estimating  error  is  to  use  energy  method        By 
equating  the  kinetic  energy  of  the  entire  system  before  impact  *o  the 
potential  energy  stored  in  the  cushioning  at  its  maximum        ".  .■■     on  the 
maximum  possible  value  of  £*J^t    can  be  determined      Tl  Is  n  -.ximum 
possible  value  can  then  be  compared  to  the  n  ix  m  in       ilue  from  the 
analog  results.    With  no  damping  (  &  l  s  ^5  6  J    and  -  i„     r  ttio    tilsUQ 
the  comparison  is  as  follows: 
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We  know  that  with    M.  -/.(),     ^/S  /^  .    Xt  =  Xx  =  ore 

impact  and      CO^  -    ft     I  vrv^     •    Substituting  the  ei 

Eq.   (8)  gives: 


,X0         -    \_   cu*      vwx    Xa  ^ 


I  »     »  i  •  i  a.  xi 


**a  x.*    =  4-  d>i  yyv,  x* 


Rearranging  gives: 


A.  = 


x 


X  /na   unM.-/ 


JL, 


Ju      (yy*^j4- 


X.  0  /  Cci 


JU 


and 


X 


«1  rw-C-J- 


Xc  /  W 


j.£ 


whi<  " 
the  maximum  possible  value  for     -jp- --:----  for  tjy,s  case.    W.^en  we 

compare  this  value  of     gi       to  the  maximum  analog  result  we  note 
that  the  analog  result  of  1.45  is  about  2.6%  high  .    7   i     .         1  error  is 
probably  a  little  higher  because  some  potential  energy  may  be  stored 
in  the  spring     }( 

A  third  method  of  estimating  error  is  to  compare  the  anal      i  al 
and  analog  results  for  the  maximum  deflections  of  the  c  .-   •      Lth 

damping  and  zero  mass  ratio.    The  maximum  deflection  of  the  cu         ning 
occurs  before  rebound  and  is  dependent  on    p      as  shown  in  Eq.   (22 
Comparison  of  the  analytical  and  analog  results  shown  in  Fi.g.  18 
suggests  of  an  accuracy  of      ±      2%„ 
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Of  further  interest  is  a  comparison  with  the  analo  j  n    lb.11     pub- 
lished in  "Dynamics  of  Package  Cushioning"  by  F     D .    ;  I  [        ho 
uses  an  amplification  factor  which  is  practic  il]     i     isl  to  tl  e    m  p 
tion  factor  used  in  this  study. 
For  perfect  rebound  with: 
/*,    =       0.01   , 

A*s       0.0    ,  .       - 

JUL    -      o.o  , 


Frequency  ratio 


1.4 

2 

2.5 

3 

4 

5 

6 

7 


Amplification  Factor 


') 


Mindlin 


1.54 
1.65 
1.61 
1.5 

1.42 
1.24 
1.08 
1.15 
1.14 
1.12 
1.06 


o& 


69 


57 


45 


07 


14 


1.13 


1,04 


The  results  compare  favorably  except  at  the  low  Ere  [uenc 
of  1.4  to  2.5  where  the  difference  is  approximately  4  -  5%. 
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Another  comparison  of  interest  can  be  made  with  analog  results  for 

inelastic  impact  published  in  "Vibration  and  Shock  Isolation     by  Co  E. 

Crede    [2]  .    Crede  reduced  data  originally  obtained  by  R.  D=  Mindlin    |Yj 

At  high  frequency  ratios  the  maximum  acceleration  of  m    occurs  before  re- 

i 

bound  so  it  is  reasonable  to  expect  the  results  for  inelastic  impact  to  be 
in  close  agreement  with  the  results  for  perfect  rebound.  For  JU^OoO  ; 
fi{  =      o.oi 


Frequency  ratio 


Trans  mi  ssibility 
Crede  Dauber 


Crede' s  values  are  taken  from  a  plot  and  are  not  to  be  considered  the 
exact  experimental  values.    It  can  be  seen  that  the  agreement  is  reasonably 
close. 

The  above  comparisons  of  analytical  and  other  analog  results  with 
the  analog  results  of  this  study  suggest  that  the  accuracy  of  the  latter  is 
approximately      t     4% . 
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Fig.   14.  Comparison  of  analytical  solutions  with  analog  results 
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Fig.  15.         Comparison  of  analytical  solutions  with  analog  results 
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Fig.  16.         Comparison  of  analytical  solutions  with  analog  results 
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Fig.   17.         Comparison  of  analytical  solutions  with  analog  results 
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V.         Conclusions 

Design  curves  have  been  obtained  for  predic    Lng  ;  of 

a  two  degree-  of  -freedom  linear  system  w  Lth  pei  nd         he  e 

curves  are  believed  accurate  within  1  4%0    The  foil  fe      ires 

of  the  response  behavior  may  be  deduced:: 

a)  Transmissibility  curves  show  that  in  order  to  I  &  efit  signific  intly 
from  cushioning  the  frequency  ratio  be    jreatei  t 

2      (Figs,  6  and  7) .    In  this  region  mere  is  in y  decreases 

.    The  addition  of  small  amount  s  of  d  m  pin  . 
in  the  cushioning  has  negligible  effect  and  jreatei  d  u  ping  increases 

(Fig.  6).    For  a  number  of  element-  g:  :    e     ame  fre- 
quency        4  ,    the  acceleration  will  be  ie  a  it  for  the  most 
massive  (largest  value  of  /*    )  (Fig.  7). 

b)  The  maximum  cushion  deflection  ie     In  ersely  with 
COx.       and  decreases  monotonically  with  incre      In  ,  damping 

(Fig.  12). 


APPENDIX  A  Summary  of  Experimental  Data 

The  summary  of  experimental  data  is  presented  in  Tables  1  to 
12. 

Tables  1  to  6  list  the  maximum  responses  with  damping  and  zero 
mass  ratio. 

Tables  7  to  12  list  the  maximum  responses  with  no  damping  and 
various  mass  ratios.    Tables  11  and  12  are  for  the  same  mass  ratio 
M  =   l.o  and  they  show  the  uncertain  scatter  established  during 
separate  runs . 
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Table  1 


0    , 

.02 

1 

'■ 

19 

■ 

.04 

38,5 

ns 

.06 

54 

constant 

&J  X    X  3  = 

.08 

68 

at  all 

.1 

78 

frequency 

!v 

.12 

82„5 

ral 

.  14  * 

86 

/       ' 

.16 
.18 

85.5 
85.5 

.2 

84„5 

/ 

.25 

.3 

.35 

78.5 
72,5 
66.5 

M 

.4 
.45 

61.5 
57 

.5 

53 

.55 

56 

.6 

57 

.65 

57 

.7 

56,5 

.75 

56 

.8 

54 

.85 

5,5 

.9 

52 

*» 

.95 

53 

1.0 

53.5 

29 


Table  2 


CO, 
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.02 

17.8 

48, 

.04 

35 

rerr  ain  s 

c 

CJL'X  -^o 

.06 
.08 

50 
62 
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at  all 

.1 

71 

frequency 

0^=0.1 

.12 

.14 

76.5 
79 

ratios 

.16 

80.5 

£  =  o-o' 

.18 

.2 

.25 

80 

78.5 

74 

! 

£*=  °'^ 

.3 
.35 

.4 

69 

64 
59 

► 

.45 

54.5 

^(    -O.Q 

.5 
.55 
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52 

.6 

53.5 

.65 

54 

.7 

54 
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53 

.8 

52 

.85 
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.95 
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51 

30 
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.02 

16 

43 

.04 

33 

remains 
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48 

constant 

U)J>.*SO 

.08 
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73 
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78 
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69 
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57 
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.5 
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52 

.6 

54 
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55 

.7 

55 
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.8 
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54 

.9 

54 

.95 

52.5 

1.0 

52 

31 


Table  4 


LO  i 

(^su    Xa 

.02 

15,2 

3/   e 

.04 

30 

remains 

o 

.06 

44 

con  itant 

&>,X.  - 

.08 

55 

at  all 

.As.     ° 

.1 

63 

frequency 

.12 

69 
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.14 

72 

<X^~-  O.I 

.16 

74 

1 

.18 

75 

.2 

75 

£,  =.  o.oi 
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74 

.3 

72 

.35 

69 

£*.=   OA 

.4 

67 

i/     **^* 

.45 

64.5 

; 

.5 

63 

=  6-o 

.55 
.6 

61 
59 

.65 

58 

,7 

59 

.75 

59.5 

.8 

60 

j 

.85 

60 

.9 

59,5 

.95 

59.5 

■ 

1.0 

59.5 

| 

Table  5 


tJ, 

CO 

x, 

\a 

,02 

13 

33 

,04 

27,5 

re:r::jr.  s 

c 

,06 
,08 

40 
50,5 

constant 

at  all 

.1 

59 

frequency 

(0*.*  oA 

,12 

65 

ratios 

.14 

69 

ft    ^o.oi 

i     i 

.16 
.18 

72 

74,5 

.2 

75 

fa-0'* 
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.3 

76 
76 

.35 

75 

Ik   =    D<° 

.4 
.45 

74 
73 

.5 
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.55 

71 

.6 

70,5 

,65 

70 

.7 

69 
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68 

.8 

67.5 
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67 

.9 
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66 

1.0 

66 

! 
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Table  6 


CJ, 

©6 

k'*    X*  m. 

.02 

13.5 

30.5 

0 

.04 

Z/  o  b 

remains 
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39.5 

constant 

.08 

50 
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.1 
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65 
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69 
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.18 

73 
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.2 
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.3 

81 
82 
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83 

JLA.  =0'° 

.4 
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.5 

84 
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83.5 

.6 

83.5 

.65 
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.7 
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Table  7 
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.7 
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.8 
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.9 
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19 
38 

53 
69 
78 
85 


88 
86 

80 

74 

68 

62 

57, 

55 

58 

59 

58 

58 

57 

56 

55 

54 

55 

55 
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at  all 
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Table  8 


U), 

NX 

u)^  XA 

",  iX 

.02 

21,5 

51  . 

u*X*Si 

.04 
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40 
56 

'  1      r 

d  JL  •  O 

.08 

68 

52 

.1 

76.5 

52 

CO^O'I 

.12 

81 

52 

.14 

83 

52 

^  =  6.0 

.16 
.18 

82.5 
82 

52 
52.5 

.2 

80 

52.5 

A  =  o.o 

.25 
.3 

74 
68 

53 

.35 

61.5 

54 

^  =  o.l 

.4 

56 

54 

.45 

53 

54 

.5 

56 

54 

.55 

57 

54 

.6 

56.5 

54 

.65 

56 

54 

.7 

55 

54 

.75 

53 

54 

,8 

52 

54 
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53 

54 

.9 
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54 

.95 

1.0 

53 
53 

54 

54 
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Table  9 


to, 

.02 

20.5 

51.5 

.04 

38 

51.5 

okX.*so 

.06 
.08 

54 
66 

\1 
52 

.1 

72 

52 

.12 

76 

52.5 

Cx).^  0.\ 

.14 
.16 

77 
76.5 

53 
53.5 

.18 

76 

54 

^o.o 

.2 
.25 
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67 

54.5 
55 

.3 

61 

55.5 
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.35 

.4 

55 
50 

56 
56 

.45 

54 

56 
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.5 
.55 

54.5 
54 
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.6 

53.5 

56 

.65 

52 

56 

.7 

50.5 

56 

.75 

50 

56 

.8 

51 

56 

.85 

51 

56 

.9 

51 

56 

.95 

50 

56 

1.0 

50 

56 
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Table 

10 
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18.4 

51 
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.06 
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50 
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.08 

60 
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.1 

65 

53 

.12 

67 
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.14 

68 

54 

£,-o.c 

.16 

66 

55 

.18 

65 

56 
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.2 

62 

58 

.25 

56 
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.3 

50 

SO 
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.35 
.4 
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49.5 
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60.5 

.45 

50 
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.6 

47 
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60 
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X 
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52 

c 
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.08 

48 

54 

.1 

50 

57 
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50 

6  0 
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.14 

48 

65 

.16 

46 

68 

.18 

44 

70 
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.2 

42 

71 
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40.5 
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.3 

42 

72 

.35 

41 

72 

.4 
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72 
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.45 

39 
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.5 
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39 

72.5 

.6 
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73 
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.7 
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38 
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38 
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47 

54 
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.12 
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49 

56 

60 

.14 

47.5 

64 
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.2 
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41 
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72 
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72 
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72 
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.65 
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72 
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.7 

38 

72 
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.75 

38 

72 

.8 
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72 

.85 
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.9 

37.5 

72 

.95 
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72 

1.0 
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| 
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Appendix  B  Analytical  Solution  to  the  Differential  Equations  oi 

Motion 
1.        Analytical  solution  for    ll<&(.0    and  no  damping  „ 
With  no  damping  Eqs .  (4)  and  (5)  become! 

X,  =   -u:/-  (  x,  -Xj .  m 

X  a.  -      -  £0  x      X  a.     —   ^  X  i    , 


The  limiting  form  of  Eq.   (10)  when  ju.    approaches  zero  is: 

X.,  =   -cJ*X 


The  initial  conditions  at  't-o    are 


x, -  x,  ,  o     . 


X,=  X,   -  X 


O       o 


Substituting  these  initial  conditions  in  the  general  solutions  of  Eqs 
(9)  and  (11)  gives: 


CO^-tJa.*    \        . 


,)<  *-  a^  L  -  u3a,   4^^  (J j  t  J  4 


Xs~  AM>f0-  Ycu,  ./Ll—a)lt'  -  c^  ^  u)^ 


U)f-tix* 


X*  =    ~£J^  AVv^  U)^t  .  (14, 

Eqs.  (12)-  (14)   are  valid  up  to  the  time  of  rebound  of  m«.    Rebound 

occurs  at  time  t    ,  when  the  acceleration  of  m^  is  zero. 

1  * 
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tlto  =      —      •  (15) 

XA--o  w^ 

After  rebound:    Losing  contact  with  the  support  requires  that  we 
replace  Eq.   (11)  with  : 

of 

XA  =    O  (ir) 

The  initial  conditions  are  determined  at  the  time  of  rebound,  "t,  ="y^' 
We  will  start  a  new  time  coordinate  with  X  -  O   .    The  initial  conditions 


are: 


^-c    w*-".*  I  Tor  u,. 


v  o 


Substituting  these  initial  conditions  in  the  general  solutions  of  Eqs. 
(9)  and  (ll1)  gives: 
X\   =■  -^Xfca>j,       <<Cfr<bTTu)t     £Vy^(  QQ^-HTU),]-  v    -£'  (16) 


J^L^L  yOHLTTfa.     ^((0^+1^).         (17) 


'3L_ 
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xA  =  -*0"t'.  <l8> 

Eqs.  (16)  -  (18)  are  valid  for    1    >  o     and    *t  ^  rt/a 
2.    Analytical  solution  for    li«\     with  damping  0 

The  differential  equations  of  motion  are  equations  (4)  and  (5) 
of  the  Introduction. 


The  limiting  form  of  Eq„  (5)  when   Ll      approaches  zero  is 
The  initial  conditions  at  u-@  are: 


,    ~  xx    -    X 


Substituting  these  initial  conditions  in  the  general  solutions  of  Eqi 
(4)  and  (19)  gives: 

x1=  eft(t  ft)+e*r(pM  ^t)t      (20) 
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+  (-  -  . •  i     -  -i  ->•  &  +  **  Ka)  e"**  a<:^  »* 

+  (  -  +-  a*  Q     -  b*  c* )  e  -**  .< ,  i   bt 

+  ( -  b^P  +  a  cbQ   +  a*  ?)  e  "•*  AU  bt . 


XA  -    X.c  e  axtw  b"t  * 


(21) 


(22) 


where  (X.  -  jo^  &>.x   .  (23  t) 


f  =   yS,  £Jy     .  (23c) 

Kl=  ^^*A^a~bp  -  (23e) 
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b  b 


Eqs.  (20)  -  (22)  are  valid  up  to  the  time  of  rebound  which  occurs 
at  time  "t,    when  the  acceleration  of  m^  is  zero,, 

-t,       =  JL  [lL  -  -u^,  -'  I  a,*-  b*  )]  .  (24) 

After  rebound:    Losing  contact  with  the  support  requires  that  we 
replace  Eq.  (19)  with: 

X^  o  (ir) 

The  initial  conditions  are  determined  at  time  of  rebound  X  fi     „    We 
will  start  a  new  time  coordinate  with  "t  -  o  <>    The  initial  conditions 
will  be  designated  as  follows: 

Substituting  these  initial  conditions  in  the  general  solution  of  Eqs  „ 
(4)  and  (19')  gives: 
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[25] 


where 

Hi  =    ^'t,  "   **•*,  •  (28a) 

Ha.=  jx,ti-  \,J  +  4(x,t|- X^() 

Eqs.   (25)  -  (27)  are  valid  for  "t'^o    and  t^t,    . 
3.        Analytical  solution  for  ft    <=.  a    _.  0       /L^o 
With  no  damping  Eqs.   (4)  and  (5)  become: 

£,    =     -  Cof*   (   X,  -    Xa-V 
Aa.  -     —  COa-    Xa  ~ycx-  Xi  « 
The  initial  conditions  at    "t  =  o      are: 

xv  =   xa *  °  j 

(a)  Because  the  solution  of  Eqs.  (9)  and  (10)  becomes  intractable  in 
literal  form,  the  following  numerical  values  are  assigned: 

Substituting  the  initial  conditions  and  the  above  assigned  values  in  the 
general  solutions  of  equations  (9)  and  (10)  gives: 
X,  =    --^AS".*!  aw.  j.^hL't  •+•  (oLS-sr,  3-,  aL^  .MHH4~t .         (23) 

X,  =    ju*  x  .^u> .  Axiu  t:  -  s"3 .  x  .iw .  ^wV  £<        (30) 
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X^=   job.s;  Auv-.Ji3.3tt  -4-5";-5o.-f  j  /u*~ .o*944t *         (31) 

Equations  (29)  -  (31)  are  valid  up  to  time  of  rebound  of    m*      .    Re-= 

bound  occurs  at  time  "t,   when  X^  is  zero.    Time  "t,    Is    letermined  by 

a  trial  and  error  procedure  and 

t:,  =^s.35"  seconds.  (32) 

After  perfect  rebound:    Losing  contact  with  the  support  requires 

that  we  replace  Eq.   (10)  with: 

f 

The  initial  conditions  are  determined  at  time  of  rebound   t-4    .    We  will 
start  a  new  time  coordinate  with    t'-c    .    The  initial  conditions  are: 
X{  ,       =  -^sr<?  x,  =  o       . 

Substituting  these  initial  conditions  in  the  general  solutions  of  Eqs . 

(9)  and  (101)  gives: 

X,r   -1$.Z-f?t..t  ¥  -+*osU^.3J9/  t'  +  4$6  Kt^.ZJf/t'.     (33) 

Kjl-  -93. SL  -*ffi.  /£'  -^ftA*^  •*•''?/'£    +73<2>  c**,2JWZ*     (35) 

Eqs.  (33)  -  (35)  are  valid  for  "£'^o    and      "£  ^      35.35  seconds, 

(b)  For  additional  comparison  with  analog  results  the  mass  ratio  jll 

is  set  at  unity  and  the  assigned  values  for    co}      uj^      and     Xd 

remain  0.2,  0.1  and  500  respectively.    Substituting  the  initial  conditions 

into  the  general  solutions  to  (9)  and  (10)  gives: 
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X,  -   -<?lS  -m^-.:-  +772  .iter.  ( 


X&-  - .  //£,*  ....  hit    *-  &g"/  8) 

Eqs.  (36  -  (38)  are  valid  up  to  time  of  rebound  which  occurs  at  time  "t 
when    X^    is  zero. 

•fc ,  =  -^6 .  o  seconds .  (39) 

After  perfect  rebound:    Losing  contact  with  the  support  requires 
that  we  replace  Eq.  (10)  with  Eq.  (10').    The  initial  conditions  are 
determined  at  time  of  rebound      "t  0    .    We  will  start  a  new  tine 
coordinate  with    t  a  o    „    The  initial  conditions  are: 


Xi  ,       =  -ns,r      t         X.   ,      -  o 

Substituting  the  initial  conditions  into  the  general  solutions  of  Eqs 
(9)  and  (10')  gives: 
X,  =  -*7«7'-T,,*W*fr# 

X,  -    /^v3     6U  .     .  +     '7-  ' 

X^=    -&7.7  -^?t       t'+/7?J*i*u.Sto        ?'+*&#  '(42) 

Eqs.  (40)  -  (42)  are  valid  for  tsK©   and    €  s  05  sec  m  I 
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Appendix  C  Summary  of  Analytical  Data 

Table  13 
Analytical  data  for: 


/ 


The  solutions  are    (  from  Eqs.  (20)  -  (22)  and  (25)  -  (27^   )    : 
Before  rebound 


y^=  ^-id*  e-^^,6„i 


After  rebound 


vct  i^  -,ooau."<r'  / 


S-5T.7      g*   -<°0*-X'   ^^  ^^ 


X.,    -     ll^o-   -~~x.^sr\?   "£" 
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Table  13 


t 

i 

t 

x, 

0 

0 

0 

3 

i 

-4.4 

13.9 

6 

-18.0 

5.1 

9 

-38.2 

32.9 

12 

-59.0 

37.0 

15 

-72.9 

37. 

18 

-78.3 

34.9 

21 

-59.2 

29.6 

24 

-31.1 

22.4 

27 

+  3.1 

14.2 

27.95 

0 

+  14.2 

11.4 

3 

+43.0 

2.8 

6 

+56.4 

-  $ o7 

9 

+50.1 

• 

12 

+26.5 

15 

-  6.0 

18 

-36.1 

21 

-53.4 

24 

-51.7 

27 

-32.3 
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Table  14 
Analytical  data  for: 


.0  / 


o.x.  *^    -    ^° 


The  solutions  are  Eqs.  (30),  (31),  (34)  and  (35): 
Before  rebound 

X, =  .  "  .  .  /  /   '  & 

After  rebound 
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Table  14 


t 

i 
t 

X, 

<*>**  x*. 

0 

0 

0 

4 

-  2.0 

19.5 

8 

-14.1 

36.0 

12 

-37.4 

47.3 

16 

-61.7 

52.3 

20 

-72.7 

51.0 

24 

-61.6 

43.8 

28 

-32.3 

31.7 

32 

35.35 

0 
2 
4 
6 

+  1.6 
+22.4 
+  29.0 
+  30.2 
+25.6 

15.5 

-  .5 

-  9.7 
-19.5 

-29.6 

8 

+  16.3 

-40.0 

10 
12 

+  3.8 

-  9.4 

-50.4 
-60.9 
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Table  15 
Analytical  data  for: 

The  solutions  are  Eqs.  (37),   (38),  (41)  and  (42): 
Before  rebound 

%\   -     ^.S'  ,aU^  ^9^/T  -  3£s,XJ<a^t0&ZS~tr4  (37) 

After  rebound 
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Table  15 


t 

i —  — 

i 
t 

x» 

"    I 

0 

0 

0 

4 

-  2.0 

19.5 

8 

-12.8 

36.3 

12 

-29.6 

49.6 

16 

-40.7 

59.5 

20 

-39.1 

66.3 

24 

-30.5 

68„8 

28 

-26.0 

65.2 

32 

-28.8 

55.5 

36 

-30.2 

41.6 

40 

-20.7 

25.8 

44 

-  2.2 

9.0 

46.05 

0 

+  7.0 

-   .  6 

2 

+  13.6 

-  9.1 

4 

+15.9 

-18.8 

6 

+  13.3 

=  29.0 

8 

+  6.6 

-39.8 

10 

-  2.3 

-50.8 
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Appendix  D  Operating  Routine 

Preliminary  preparations:    The  coefficient  potentiometers  were 
calibrated  with  the  comparator  bridge.    This  calibration  was  verified 
several  times  during  the  course  of  the  study.    The  recorder0     timer  was 
calibrated  with  an  electric  stop  watch. 

Daily  preparations:    The  computer's  amplifiers  were  balanced 
after  an  initial  warmup  period  of  at  least  twenty  minutes  and  at  the 
beginning  of  each  series  of  runs.    The  initial  conditions  were  verified 
after  each  series  of  runs .    The  recorder  channels  were  calibrated  before 
each  series  of  runs.    The  calibration  was  checked  after  the  series  of 
runs  and  recalibrated  if  necessary.    Capacitors  and  resistors  were  cali- 
brated with  a  capacitance  bridge  and  a  standard  resistor  respectively. 
Each  capacitor  was  inserted  in  a  separate  integrating  circuit  and  checked 
for  improper  period  and  excess  leakage.    Coefficient  potentiometers  in 
the  problem  circuit  were  adjusted  for  any  day  to  day  variation  of  the 
capacitors.    A  check  off  sheet  was  maintained  for  all  daily  preparations. 

In  all  runs  the    CUg     potentiometers  were  kept  at  a  constant  value 
of  0.1  rad.  /sec.  and  the   &>,    potentiometers  were  varied  to  give  the 
desired     CJf  / CU^    ratio.    The  ftt     ^    and  fi    potentiometers  were  set 
with  their  required  values  and  remained  constant  for  the  complete  range 
of  the      Oit  /  ratio.    The  output  of  Amplifier  #1  was  always  fed  to 

one  of  the  recorder  channels  to  verify  proper  operation  of  the  diode 
limiter.  At  cutoff  by  the  diode  limiter  the  output  voltage  of  Amplifier  #1 
was  approximately  0.4  volts.    The  maximum  output  of  Amplifiers  #5  and 
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#6  were  recorded  on  the  data  sheets  „    A  rough  plot  of  ftr 

versus    UJjju! ,   was  made  concurrent  with  computer  operation  „    This 
rough  plot  served  as  an  indication  of  unexpected  responses. 
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Appendix  E  Equipment 

a.  Analog  Computer,  Model  30 
Donner  Scientific  Company 
Berkeley,  California 

b.  Comparator  Bridge  ,  Model  50 
Donner  Scientific  Company 
Berkeley,  California 

c.  Four  Channel  Recorder,  Model  154-1 100A 
Sanborn  Company 

Cambridge,  Mass. 

d.  Vacuum  Tube  Voltmeter ,  Model  311 
Simpson  Electric  Company 
Chicago,.  Illinois 
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